Minimal self-contained quadridot quantum refrigeration machine 
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We present a theoretical study of an electronic quantum refrigerator based on four quantum dots 
arranged in a square configuration, in contact with as many thermal reservoirs. We show that the 
system implements the basic minimal mechanism for acting as a self-contained quantum refrigerator, 
by demonstrating heat extraction from the coldest reservoir and the cooling of the nearby quantum- 
dot. 
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The increasing interest in quantum thermal machines 
has its roots in the need to understand the deep relations 
between thermodynamics and quantum mechanics 0, Q ■ 
The progress in this field may as well have important 
applications in the control of heat transport in nano- 
devices Q . In a series of very interesting works the 
fundamental limits to the dimensions of a quantum re- 
frigerator have been found. It has further demonstrated 
that these few-level quantum machines could still attain 
Carnot-efficiency [B| thus launching the call to find phys- 
ical implementations of the smallest possible quantum 
refrigerator. Refs.jH-Q considered self-contained thermal 
machines defined as those thermal machines that per- 
form a cycle without the supply of external work. Their 
action being grounded on the steady-state heat transfer 
from thermal reservoirs maintained at different temper- 
atures. While such devices are ordinary in the macro- 
scopic world, their realization in the quantum realm is 
challenging. 

As recognized in recent works 0, Q on self-contained 
refrigerators (SCR) the major difficulty in the realization 
of the three-qubits quantum machines of 0, [H[ is the engi- 
neering of the crucial three-body interaction enabling the 
clastic transition between a doubly excited state in con- 
tact with a hot (H) and cold (C) reservoir, and a singly- 
excited state coupled to an intermediate-temperature (or 
"room-temperature" - R) bath. In this paper we get 
around this problem by proposing an experimentally fea- 
sible implementation of a minimal quantum refrigerator 
in a small array of semiconducting quantum dots (see 
Fig. Q]) operating in the Coulomb blockade regime. With 
this work we are thus able to establish a further con- 
nection between the general theory of quantum thermal 
machines and the heat transport in nanoelectronics 0]. 

QDs contacted by leads have been recognized as proto- 
typical systems for achieving high thermopower, see 0- 
ll| and references therein, or anomalous thermal ef- 
fects [13] ■ Here we study a four-QD planar array (here- 
after named a "quadridot" for simplicity) coupled to four 
independent electron reservoirs as shown in Fig. [TJ with 
proper (but realistic) tuning of the parameters, we will 
show that this system can effectively represent an implc- 




FIG. 1: [Color online] The quadridot. Vertical and Horizontal 
interdot capacitances are considered to be the highest energy 
scales among those in the figure. Their values U\\ and U± 
are determined by the arrangement of the top gates over the 
QDs, which are not shown. The four quantum dots QCh, 
QC C , QCr 1 are QCr 2 are weakly tunnel coupled to the elec- 
tronic reservoirs H, C, Ri, and R2, respectively, which are 
all grounded but maintained at equilibrium at a well-defined 
temperatures Th > (Tr 1 — Tr 2 — Tr) > Tc- No tunneling is 
possible between dot 1 (orange) and dot 2 (blue), and between 
dot 4 (red) and dot 3 (orange). [See text for the definition of 
the other symbols]. 



mentation of the model described in Ref. |j] . Specifically 
we will show that the quadridot acts as a refrigerator 
which pumps energy from the high temperature reservoir 
H and the low temperature reservoir C to the intermedi- 
ate temperature reservoirs Ri , R2. Furthermore we will 
analyze the conditions under which the quadridot is able 
to cool the quantum dot QD C which is directly connected 
to the bath C, at an effective temperature that is lower 
than the one it will have in the absence of the other 
reservoirs. This will lead us to introduce an operative 
definition of the local effective temperature depending 
on the measurement setup, and to predict the existence 
of intermediate working regimes where, for instance, the 
refrigeration is not accompanied by the cooling of QD C . 
In deriving our results we start from a detailed analy- 
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FIG. 2: [Color online] Pictorial view of the low-energy elec- 
tronic charged states of the quadridot of Fig. [Tj when it is 
isolated from the thermal reservoirs (a black circle inside the 
dot indicate that the latter is occupied by an electron. The 
completely empty state of energy E° — is not shown). The 
indicated energies correspond to the t = case. The two 
electron states \d) and \d) are resonant, while \u) and \l) are 
the high-energy virtual states responsible for the effective cou- 
pling g dl i which couples them after the unitary transformation 
(see text). The remaining 7 high-energy states (2 two-particle 
states, 4 three-particle states and 1 fully occupied state) are 
not drawn as their contribution to the projected dynamics is 
negligible. 



sis of the system Hamiltonian to identify the working 
regimes that allows us to mimics the behavior of the SCR 
of Ref. Q. 

In the absence of the coupling to the leads, the array 
of QDs (the quadridot) shown in Fig. [T]is described by 
the Hamitonian 

T~Lqd = e i n i +'^Uijninj-t(clc4+clc3+h..c.), (1) 

1=1.4 ift 

where ei,... ,4 denote the single particle energy level of the 
dots, and where we introduced the creation/annihilation 
operators cj, Cj of an electron in the i-th dot and m = 
cjcj. The matrix Uij takes into account the finite-range 
contribution of the Coulomb repulsion while the parame- 
ter t gauges the tunneling coupling. The level separation 
in each dot is assumed to be large enough such that it 
is sufficient to consider only one level. Moreover the on- 
site repulsion is taken as the largest energy scale in the 
problem reducing the maximum occupancy in each to 
to one electron. Gates are applied to the dots in order 
to tune their single-electron energy level of the upper- 
right dot (which will be coupled to the cool resorvoir) to 
€2 = ei + £3 — 64. We assume that U12 = U34 = U± and 
U23 = U14 = U\\, both much larger than the diagonal 
term U24 = U13 = Ud- The charge states of our system 
are labeled according to the occupation of the four dots 
|ni, ri2, TI3, 714). The single electron level £2 is chosen so 
that in the absence of tunneling (t = 0), the "diagonal" 
two-particle states \d) = |1, 0,1,0) and \d) = |0, 1,0,1) 
shown in Fig. [2] are, for t = 0, degenerate (in our model 
they play the role of the states |01), 1 10) of Q). 

Before going in the detailed description of the refrig- 
erator, it is instructive to see what are the states that 
will play a leading role in the thermal cycle. The low- 
est energy subspace can be identified by the four states 



with single occupancy |1) = |1, 0,0,0), |2) = |0, 1,0,0), 
|3) = |0, 0,1,0), and |4) = |0, 0,0,1). In the two-electron 
sector only the states \d) and \d) are relevant. The two 
other states \u) = |1, 1,0,0) and |Z) = |0, 0,1,1) arc 
much higher in energy, as they involve U± or U\i. They 
however play a fundamental role in the refrigerator as 
they allow to generate and effective coupling between 
the states |d) and \d). This can be see by performing 
a Schricffcr- Wolff transformation [13 , 141: one can show 



that besides a slight-renormalization of the single parti- 
cles levels (which we will ignore) , a direct coupling among 
\d) and \d) emerges with an intensity which can be well- 
approximated by 



9dd 



2t 2 (U d -U ± ) 



(U d - U±_f - (e 4 - ei) 2 



< t. 



(2) 



In our model this is the analogous of the crucial pertur- 
bative g term discussed in Ref. 0] as it allows to per- 
form a coherent swap between the state connected to 
the room-temperature reservoir (\d)) to the doubly oc- 
cupied state \d) which is simultaneously in contact with 
the hot and cold bath. Due to the effective hoppings of 
t and g d d the eigenstates of the low-energy Hamiltonian 
will be the bonding-antibonding states \d) ± \d) and the 
four bonding-antibonding single particles states delocal- 
ized between the dots 1, 4 and 2, 3. 

As anticipated, each QD is contacted to an indepen- 
dent electron reservoir characterized by a chemical po- 
tential fj,i and temperature Ti, appearing in the Fcrmi- 
Dirac occupation functions /j(e). For the sake of the 
simplest correspondence with the model of Ref. 0] in 
this study all [ii will be set to be identical and fixed 
to a value that will be used as reference for the sin- 
gle particle energies of the system [e.g. setting = 
in Eq. ([1]) corresponds to have the i-th quadridot 
level at resonance with the Fermi energy of the reser- 
voirs]. Furthermore the temperatures will be chosen as 
T, = T 3 = T R , T c = T 2 and T H = Ik satisfy T c < 
Tr < Th for the refrigerator to work [15| [the change 
in the notation helps in identifying the cool, hot and 
room temperature reservoirs]. Within these assumptions, 
the only required external action is exerted in main- 
taining the local equilibrium temperatures and chemi- 
cal potential, in accordance to the standard definition of 
self-containancc. The equilibrium correlation functions 
(ai.k/a\ fe / destroys/creates an electron in lead i and mo- 
mentum k) read J^k I e%et ( < A,k( t ) a i,k{Q))dt = / i (e)^(e) 
and EkI e l£t (a hk (t)al k (0))dt = (1 - /i(e))^(e). In the 
previous expression v(e) is the local density of states of 
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The coupling to the leads (parametrized 



by the amplitudes F^) reads 



££lf>cta iifc + h.c. 

i=l,4 k 



(3) 



giving rise, in the Born-Markov-Secular limit 171 [l_8(, to 
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a Lindblad master equation for the reduced density ma- 
trix p of the quadridot. The presence of t plus the unitary 
rotation into the low-energy sector eliminates any degen- 
eracy between all possible energy transitions between the 
eigenstates of the quadridot. The evolution of p is then 
determined by the master equation 



P 



(4) 



where for each reservoir T>i represents associated Lind- 
blad dissipator. This equation can be solved in the 
steady-state regime (p = 0) yielding the asymptotic con- 
figuration p°°, from which the heat currents (JQ,i) flow- 
ing through the i-th reservoir are then easily computed 
as 0, 



Tr Hqd Vi 



(5) 



If the correspondence between our device and the theo- 
retical qubit model of [H is correct, we should see a direct 
heat flow from the hot and cold reservoirs to the room 
temperature reservoirs, while dot coupled to the cold 
reservoir should reach an occupation probability lower, 
£2 > than the one dictated by its local reservoir (See 
Fig.^b). We have verified this by setting the system pa- 
rameters to be consistent with those presented in Ref. [H 
(making sure however that for such choice no additional 
degeneracies are introduced into the system due to the 
larger dimension of our physical model). For this reason 
we choose e a = 2.1, e^= 2.9, e 4 = 4.0, T, = 0.001 (mo- 
mentum independent [la ]. U± = 12.0, while kept U\\ = oo 
for simplicity as its effect could be absorbed in the energy 
level renormalization after the Schrieffer- Wolff transfor- 
mation. The value of g d( j is then fixed to be -0.001 deter- 
mining t (< 0.1) through Eq. ^ (in all these expressions 
the energy units should be considered as multiples of 100 
peV in order to be experimentally reasonable). Solving 
numerically the steady state equation (j4]) we then notice 
that for each Tc < Tr, there exists a minimal threshold 
for Tjj above which the refrigerator indeed extracts 
heat from the cold reservoir C. The results are plotted in 
Fig.[3]-a for Tr=2 and different values of U d : notice that 
consistently with the second principle of thermodynam- 
ics Tjj is always greater than Tr = 2 (for Tjj below Tr 
the machine cannot produce work from H to pump heat 
from C), and that the functioning region Tjj > Tjj gets 
larger as U d gets smaller. The existence of the thresh- 
old Tjj implies also that, for given Tjj > Tr, there is 
a minimal temperature T c for the cold reservoir under 
which the refrigerator cannot work. Interestingly enough 
for Tjj /Tr — > oo the value of T c appears to asymptoti- 
cally converge toward a finite non-zero temperature T c 
which depends upon the engine microscopic parameters 
and which can be interpreted as the emergent absolute 
zero of the model. An approximate analytical expres- 
sion for T c (and hence Tjj) can be derived exploiting the 




U d =2.0 



0.0 0.5 1.0 1.5 Tc O.O 0.5 1.0 1.5 Tc 0.0 0.5 1.0 1.5 Tc 

b > * _U d . z . 4) 

TiiT 

c 




0,0 0,5 LD 1.5 Tc 2 4 8 8 10 12 Th 



FIG. 3: [Color online] a) Working conditions of the refriger- 
ator. Panels refer to Ud — (left), Ud = 1 (center), Ud = 2 
(right). The change of sign in the heat flow occurs at T^ at 
a value very close to the analytical value of Eq. ([6]) (dashed 
line). For Th > T^ (blue region) machine works as pictured 
in b), extracting heat from the C,H-reservoirs and pumping 
it into R. Black dashed line above the grey region indicates 
Th = Tr=2. Blue/Red background color intensity is pro- 
portional to the actual heat pumped to/extracted from C- 
reservoir. c) Possible temperature regimes for Ud = 2. In re- 
gions I /II the refrigerator is working (heat is extracted from 
C) while in regions III/IV the C bath receives heat. In regions 
I/III we have an effective decrease of single particle occupa- 
tion number (i.e. no < n^). d) Efficiency of the refriger- 
ator compared to the Carnot efficiency (dashed black line) 
for Tc=1.0. The curves represents Ud = (blue), Ud — 1 
(magenta) and Ud — 2 (brown). The dashed horizontal line 
indicates maximum limit of efficiency for the quadridot com- 
puted (for Ud = 0) as in Q]. 



recent general theory of genuine, maximally-efficient self- 
contained quantum thermal machines This is done 
by interpreting the quadridot as a composite system, con- 
sisting of the QD C at equilibrium with its local reservoir, 
in contact (by means of g d( [) with an "effective" virtual 
qubit whose states are defined by |0, 0, 1,0) and \d) (See 
Fig. [2]) and whose average occupation is determined by 
a virtual temperature that can be interpreted as our Tq. 
Following this yields the expressions 

T^T R T H TH{ lXZ- + ^ TRei , T^T R (1- K ), (6) 

with k = + £3 + U d ], which fit pretty well our nu- 

merical results (See Fig. [5]-a) fl9j . 

Following we also evaluate the ratio r\ = 

(Jq,2)/(Jq,4) between the heat current through the cold 
and hot reservoirs comparing it with the upper bound 
(1 — — 1) posed by the Carnot limit, and with 

the theoretical value of the three-qubit model [ij applied 
to the quadridot for U d = [§oj], r)th={<^i + £3 — £4)/e4- 
The dependence of 77 as a function of Tjj arc plotted in 
Fig. El-d for different values of Ud- We noticed that in the 
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special case Ud=0 the efficiency of the quadridot con- 
verges indeed towards the theoretical value r\th of 3 at 
least for large enough T H . 

Effective Local Temperatures:- The sign of the heat 
flow in the working region suggests that any object (de- 
fined by a local Hamiltonian) interposed between the 
QD2 and its bath will necessarily be cooled by the SCR. 
An important question is whether this refrigeration effect 
is accompanied with a cooling of the interacting QD C , 
namely whether its effective local temperature Tq ' de- 
creases as Tjj increases, for sufficiently high TW, in anal- 
ogy with the qubit-cooling described in Ref. [4|. While 
for such idealized qubit model the definition of the local 
temperature is relatively straightforward, in nanoscale 
systems out of equilibrium, local temperatures must be 
operationally defined 2l| . A natural way to proceed is to 
introduce a probe reservoir P (a "thermometer") which 
is weakly coupled to that part of the system we are in- 
terested in (the dot QD C in our case) and identify the 
effective temperature of the latter with the value of the 
temperature Tp of the probe which nullifies the heat flow 
through P. This procedure can be implemented easily in 
our model by adding an extra term in that connects 
the new reservoir P to QD C with a tunnel amplitude Tp 
which is much smaller than those associated with the 
other reservoirs of the system (in the calculation we set 
the ratio between Tp and Tc,t,r to be of the order 10~ 3 : 
this make sure that the presence of P does not perturb 
the system) . The obtained values of are presented 

in Fig. [3}a where it is shown that, according to this defi- 
nition of the local temperature, the conditions for cooling 
of QD C (i.e. Tq ' < To) are the same for the refriger- 
ator to work (in particular this implies is always 

greater than the emergent zero-temperature of the sys- 
tem). 

The quantity introduced above has a clear oper- 

ational meaning and it is our best candidate to define the 
effective temperature of QD C . Still it is important to ac- 
knowledge that in experiments the cooling of QD C could 
also be detected by looking at the decrease of the mean 
asymptotic occupation number of the QD C (i.e. (nc) = 
(0, 1,0, 0|p°°|0, 1,0,0) + (d\p°°\d)), with respect to the 
same quantity computed when the refrigerator is "turned 
off" (e.g. (n° c ) = (0, 1, 0, 0|pg°|0, 1, 0, 0) + (d|pg°|d) where 
now is the asymptotic stationary state of the system 
reached when all the reservoirs but C are disconnected, 
i.e. Th = = Tp 2 = 0). Wc notice however that the 
cooling condition hereby defined does not coincide with 
the same pictured in Fig. [3}a. We indeed exemplify in 
Fig. [3}c for Ud=i that according to this new definition 
of cooling different operating regimes are generally possi- 
ble for the refrigerator. The QD C might be either colder 
(nc < n c in zone I) or hotter (nc > n c , in region 
II) when the device extract heat from the C reservoir. 
Conversely, we might achieve a colder QD C also when 



the quadridot pumps heat into the colder bath (III). For 
the temperatures region IV none of the refrigeration ef- 
fects is active. For completeness, we also report that 
the effective local temperatures of the other three QDs 
are generally positively correlated with the increase of 
Th, witnessing the intuitive consideration that cooling 
of one part of the quadridot is made at the expenses of 
the heating of the rest of the device, independently from 
the operational definition of local temperature. Similar 
intermediate regimes emerge with other activation pre- 
scriptions, such as defining (n c ) as the occupation for 
Th = Tp : — Tc while maintaining all tunnel couplings as 
constant. 

Conclusions:- We conclude by commenting on the 
possibility of experimentally attaining the physical 
regimes described in the manuscript. Quadridots in stan- 
dard GaAs/AlGaAs heterostructures have been imple- 
mented aiming at Cellular- Automata computation [22 1 
and very recently for controlled single-electron manipu- 
lation [231 ]. Strongly capacitively-coupled QDs with in- 
terdot capacitance energy (U± and Uu) up to 1/3 of the 
intra-dot charging energy (taken to be infinite in our 
model) can be fabricated with current lithographic tech- 
niques [2^ | . The diagonal inter-dot term Ud is expected 
to be at most U\\/y/2 ~ U±/\/2 from geometrical con- 
siderations, but practically it is expected to be much 
smaller [23[ and in principle can be reduced arbitrar- 
ily with appropriate dielectrics or by sputtering in the 
center of the quadridot. Experimentally, the local charg- 
ing energy can be as big as 1 meV, and usually repre- 
sents about 20% of the confinement energy [25[ , which is 
the typical tunable values of the single-particles levels . 
While these ratios can already partially justify the spec- 
tral scales necessary for the refrigerator to work, charging 
effects (relative to single-particles levels) are expected to 
be further enhanced by the presence of a significant mag- 
netic field, due to the emergence of the incompressible 
antidot regime in the dots }26j , possibly allowing the nec- 
essary working conditions to be achieved more easily. We 
also note that in this high-field regime, the spin/orbital- 
Kondo effect (which is expected active in these kind of de- 
vices [I?], [28j| ) is suppressed as the transport becomes 
spin-polarized, so our effective description is expected to 
be valid. We finally note that the maximum thermal 
energies involved should not exceed the large charging 
energies, so the bath temperatures should all be < 10K. 
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IDEAS project RBID08B3FM, by EU through the 
Project IP-SOLID, and by KK9. 
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